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Abstract 

The dynamics of Fabry-Perot cavity with suspended mirrors is described. The suspended 
mirrors are nonlinear oscillators interacting with each other through the laser circulating in the 
cavity. The degrees of freedom decouple in normal coordinates, which are the position of the 
center of mass and the length of the cavity. We introduce two parameters and study how the 
dynamics changes with respect to these parameters. The first parameter specifies how strong 
the radiation pressure is. It determines whether the cavity is multistable or not. The second 
parameter is the control parameter, which determines location of the cavity equilibrium states. 
The equilibrium state shows hysteresis if the control parameter varies within a wide range. We 
analyze stability of the equilibrium states and identify the instability region. The instability is 
explained in terms of the effective potential: the stable states correspond to local minima of the 
effective potential and unstable states correspond to local maxima. The minima of the effective 
potential defines the resonant frequencies for the oscillations of the cavity length. We find the 
frequencies, and analyze how to tune them. Multistability of the cavity with a feedback control 
system is analyzed in terms of the servo potential. The results obtained in this paper are general 
and apply to all Fabry-Perot cavities with suspended mirrors. 
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1 Introduction 



Very long Fabry-Perot cavities serve as measuring devices for interferometric gravitational wave de- 
tectors. Several such detectors are currently under construction [Jl], [3]]. The cavities are planned 
to have high circulating power and large storage time. For example, LIGO (Laser Interferometer 
Gravitational wave Observatory) Fabry-Perot cavities will accumulate 10 kW of power and will have 
roughly 1 ms of storage time. 

The suspended mirrors are designed to move freely along the direction of the beam propagation. 
Due to multi-beam interference in the Fabry-Perot cavity the motion of the mirrors strongly affects 
the light inside the cavity. The light, in turn, affects the motion of the mirrors by exerting radiation 
pressure on them. The interaction of light in the cavity with the suspended mirrors through radiation 
pressure gives rise to a nonlinear dynamics. Finite time of light propagation in the cavity introduces 
a time delay in the dynamics. The propagation time gives rise to storage time. Thus a Fabry-Perot 
cavity is a dynamical system with delay; such systems are known to have instabilities [Q]. 

The significance of the "spring" action and the "damping" effect of the radiation pressure for the 
dynamics of the Fabry-Perot cavity was pointed out by Braginskii [|5p. The cavity with one suspended 
mirror driven by radiation pressure was studied experimentally by Dorsel et al []6|, [7J, ||, §]. The 
main results reported in these papers are observations of optical bistability and mirror oscillations 
with frequencies determined by the radiation pressure. These authors analyzed their results using 
the adiabatic approximation for the intra-cavity field. At about the same time Deruelle and Tourrenc 
studied the Fabry-Perot cavity with suspended mirrors theoretically [ |10[ p| . Their analysis revealed 
delay-induced instability caused by the radiation pressure in the cavity. This instability was further 
studied by other researchers [ |I2| , [13[ ]. Stabilization of the Fabry-Perot cavity by a control system 
was discussed by Meers and MacDonald JT4|]. Recently, the radiation pressure induced dynamics 
of Fabry-Perot cavities attracted attention of the VIRGO group in connection with the design of the 
length control system of their detector [[TJJ]. Similar research has been done in LIGO and is presented 
in this paper. 

Emerging laser gravitational wave detectors require detailed modeling and pose new questions for 
the study of dynamics. From a phenomenological point of view, there is a question of what are the 
parameters that define the universal properties of Fabry-Perot cavities with suspended mirrors, and 
how the dynamics changes with respect to these parameters. From a point of view of applications, 
there is a question of how to generalize the results obtained in table-top experiments to large scale 
Fabry-Perot cavities of the gravitational wave detectors. 

In this paper we attempt to provide a phenomenology of the Fabry-Perot cavities for modeling 
and optimization of the performance of LIGO interferometers. Due to the complexity of the subject 
we split the discussion into two papers. In the first paper we study various aspects of the nonlinearity 
in the dynamics, leaving aside the time delay instabilities. In the second paper [ |T6| ] we consider the 
time delay instabilities and study the dynamics of a Fabry-Perot cavity with a realistic control system. 

In this paper we formulate the dynamics in terms of normal coordinates: the cavity length and the 
cavity center of mass. We show that a small part of the radiation pressure in the cavity, the excess 
force, excites the cavity center of mass. In absence of the excess radiation pressure, the dynamics 
of the cavity length is equivalent to the dynamics of a suspended mirror in a cavity, which has one 
mirror suspended and one mirror fixed. To study the universal properties of the cavity dynamics, 
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such as multistability, we introduce two parameters. The first parameter is a control parameter which 
allows us to select the equilibrium state. The second parameter characterizes strength of the radiation 
pressure and determines whether the system is stable or multistable and how many equilibrium states 
are there. The results obtained in this paper are general and apply to any Fabry-Perot cavity with 
suspended mirrors. Numerical calculations and modeling with parameters of LIGO cavities are given 
throughout this paper. 

The paper is organized as follows. In Section 2 we describe the equations of motion and define 
the control parameter. In Section 3 we formulate the dynamics in normal coordinates. In Section 4 we 
construct equilibrium states and introduce the specific radiation pressure. In Section 5 hysteresis and 
instabilities are described. In Section 6 the global properties of the dynamics are analyzed in terms of 
the effective potential. In Section 7 we provide explicit formulas and numerical values for the cavity 
resonant frequencies. In Section 8 we briefly describe nonlinear dynamics of the cavity with a control 
system. 
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2 Dynamical Equations and Regimes 



2.1 Mirror Coordinates 

The dynamical system we study in this paper is a Fabry-Perot cavity with two suspended mirrors, 
labeled by a and b, and a laser incident on the cavity from one side. A suspended mirror is a rigid 
body with six degrees of freedom, whose dynamics depends on the suspension design. In this paper 
we neglect angular degrees of freedom of the mirrors and consider only the motion of their center of 
mass. We also neglect the side motion of the mirrors. In this approximation the mirrors are always 
aligned to the incident beam and move parallel to themselves as simple pendula. We specify the 
mirror positions by their coordinates, x a (t) and £&(t), in the inertial coordinate frames as shown on 
figure [p. The mirror suspension points are also moving, their coordinates are x' a (t) and x' b (t). 



Figure 1 : Coordinates of mirrors and their suspension points 




reference o reference 

plane a plane b 



The two inertial coordinate frames are separated by a large distance, the nominal cavity length, 
L , which we consider fixed. The actual length of the cavity depends on mirror positions in their 
respective coordinate frames: 

L = L + x b (t) — x a (t). (1) 
The large cavity length of the gravitational wave detectors makes the delay time, 

T = — , (2) 

c 

to be non-negligible. For example, in LIGO interferometers the cavities are 4km long and the corre- 
sponding delay time is roughly 13 fj,s. 
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2.2 Equations for Field Dynamics 



Multibeam interference of laser in long Fabry-Perot cavities of gravitational wave detectors is far 
more complex than a similar process in short Fabry-Perot interferometers with fixed length which are 
commonly used for laser spectroscopy. Both the motion of suspended mirrors and the delay times 
affect the interference and give rise to a complex dynamics. We describe dynamics of the laser in the 
cavity in terms of complex amplitudes of the electric field. These complex amplitudes correspond to 
the traveling electro-magnetic waves and are defined at some locations. Usually the amplitudes are 
defined at the mirror surfaces. However, the suspended mirrors are moving and such a choice would 
require account for Doppler effects. Instead, we define the amplitudes at the reference planes which 
are fixed in the inertial coordinate frames. 

Figure 2: Notations for fields and reference planes 

reference reference 
plane a plane b 




Our notations for the amplitudes of the electro-magnetic waves are shown in Fig. ^. The amplitude 
of the wave, incident on the front mirror, is E in . The amplitudes of the forward-propagating wave 
inside the cavity, defined at the two reference planes, are E 1 and E[. Similarly, the amplitudes of the 
backward-propagating wave are E 2 and E' 2 . The amplitudes of the cavity transmitted and reflected 
waves are E 3 and E 4 . 

Propagation of the laser inside the cavity is described by the delay equations: 

E[(t) = E 1 (t-T), (3) 
E' 2 (t) = E 2 (t-T). (4) 

These and other equations for the complex amplitudes are derived in the Appendix |B]. The equations 
for reflection and transmission of the laser at the end mirror are: 

E 2 {t) = -r b e- 2ikXb(t) E[{t), (5) 
E 3 {t) = t b E[(t), (6) 

where r b is the reflectivity and t b is the transmissivity of the end mirror. 

Using the superposition principle we obtain similar equations for the front mirror: 

E^t) = t a E m -r a e 2ikx ^E' 2 (t), (7) 
E 4 {t) = t a E' 2 {t)+r a e- 2%kx ^E m , (8) 
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where r a and t a is the reflectivity and the transmissivity of the front mirror. Numerical values for the 
reflectivities and the trasmissivities of the mirrors in LIGO cavities are given in the Appendix 

The equations (|3|-|8|) describe dynamics of the laser interacting with the moving mirrors of the 
cavity. These equations can be reduced to one equation with one field, which is called the self- 
consistent field. All other fields can be expressed in terms of the self-consistent field through the field 
equations above. 

We choose the forward-propagating field to be the self-consistent field. Let the amplitude of the 
self-consistent field be E(t) (same as E\(t)). Then the equation for dynamics of the self-consistent 
field is 

E(t) = t a E m + r a r b e- 2lkz ^E{t - 2T), (9) 
where z(t) is the relative mirror position defined as 

z(t) =x b (t-T)-x a (t). (10) 

The quantity, L + z(t), is the distance between the mirrors in relativistic sense. 

Since the speed of light is finite we cannot obtain information about the instantaneous cavity 
length defined in eq. ([I]). However, we can measure the relativistic cavity length, z(t), by observation 
of one of the cavity fields. This can be done using the Pound-Drever signal extraction scheme. The 
Pound-Drever signal is based on observation of the cavity reflected laser 

Vp. D ~ lm{e 2ikXa{t) E 4 {t)} (11) 
~ Im{£(t)}, (12) 

which is defined in terms of the self-consistent field and, therefore, depends on the cavity length, z(t). 

In the expression, eq. ([K]), the delay appears in the coordinate of the end and not the front mirror. 
This is because upon entering the cavity the laser first reflects off the end mirror and then reflects off 
the front mirror. 



2.3 Forces Acting on Mirrors 

Dynamics of the suspended mirrors is defined by the forces acting on them. Since we analyze only 
longitudinal motion of the mirrors we consider only horizontal components of the forces. The hori- 
zontal component of the wire tension, acting on the front mirror, is 

-m a u%(x a - x' a ), (13) 

where m a is the mass of the front mirror and ujq is the pendulum frequency. There is a similar 
expression for the end mirror. 

The laser circulating in the cavity and the incident laser exert pressure on the mirrors. Let the 
radiation pressure on the front and the end mirror be R a and R b . These forces are defined by the 
amplitudes of the fields as follows: 

R a = - c {\E in \ 2 -\E 1 \ 2 -\E' 2 \ 2 + \E A \ 2 }, (14) 

R b = -{l^r + l^l 2 -!^! 2 }- (15) 



7 



Many control schemes, including the control system of LIGO arm cavities, require actuators on 
the mirrors. The actuators can be coil-magnet pairs, electro-static boards or other devices. For our 
analysis the details of the actuators are not important. We only assume that actuators are perfectly 
aligned and apply forces to the mirrors in the direction of the beam propagation. Let the actuator 
forces be F% ct and F b act . 

There is also a friction force due to energy losses in the suspension wires and the actuators. We 
assume that the friction is proportional to velocity. Then the friction force acting on the front mirror 
can be written as 

- m a -fi a , (16) 

where 7 is the damping coefficient. A similar expression can be obtained for the end mirror. Typically 
the damping coefficient is very small: 'j/coo ~ 1CT 5 . 

Numerical values for the parameters of LIGO interferometers are given in Appendix [A|. 



2.4 Control Parameter 

Let us consider the sum of the pendulum restoring force and the actuator force for the front mirror: 

-m a ul(x a -x' a )+F: ct . (17) 

There are two parameters here: the position of the suspension point and the actuator force. The mirror 
can be controlled through both the actuator and the suspension point. However, the dynamics depends 
only on the sum of the two parameters, and not separately on each one. Therefore, we combine the 
two parameters into a new parameter 

Ua = X a + (18) 

This is the control parameter for the front mirror. A similar parameter, tt 6 , can be introduced for the 
end mirror. 

The control parameter has dimension of length and the meaning of an equivalent displacement of 
the suspension point. The equivalence is based on the conversion factor (mw 2 ) -1 , which for LIGO is 
roughly 4.2 millimeter per newton of the actuator force. 

Moving the suspended mirrors at high frequencies is easier using the actuator. However, applying 
large displacements to the mirror at low frequencies is better through the suspension point. Therefore 
in a realistic control system the control parameter, u, may obtain its high frequency component from 
the actuator force and its low frequency component from the position of the suspension point. 

2.5 Equations of Motion for Mirrors 

The equations of motion for the mirrors are 



X a + !X a + UJ^Xa = U%U a + — , (19) 

m a 
R 

Xb + ix b + ulxb = ujnUb H . (20) 
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These are the equations for a pair of damped harmonic oscillators driven by forces of the radiation 
pressure. The oscillators are not independent: they are interacting with each other through the laser 
which acts as a spring connecting the mirrors. The "spring" is nonlinear and has time delay. 



2.6 Dynamical Regimes 

A Fabry-Perot cavity has two intrinsic time scales: one is the delay time, T, the other is the storage 
time 

2T 

T = TT- ( rr. (21) 

I m{r a r b )\ 

Correspondingly, there are two intrinsic frequencies defined by the cavity: the free spectral range 
(/fsr = (2T) _1 ) and the cavity low pass frequency 

fcav 7> ' (22) 
27TT 

Therefore, there are three dynamical regimes for the Fabry-Perot cavity. 

Let St be a characteristic time for a dynamic process involving changes in the mirror positions and 
fields. Let / be the characteristic frequencies of the process. The three regimes are: 



1. decoupled regime: St <C T (/ >> /fsr), 

2. delay regime: T < St < r (f cav < / < /fsrX 

3. adiabatic regime: St 3> r (/ 3> / ca «)- 



The first regime corresponds to very fast changes in mirror positions and fields in the Fabry-Perot 
cavity. The decoupled regime can be understood on the following example. First the mirrors are at 
rest. Then one mirror acquires some velocity, changes its position and stops. If the entire process 
takes time less than T the second mirror becomes affected by the process after the first mirror comes 
to rest. Therefore, during such process the mirrors are independent of each other. This is nothing but 
a relativistic causality limit. In LIGO interferometers the processes with frequencies above 37.5 kHz 
belong to the decoupled regime. 

The second regime takes place when the changes in the mirror position are faster than transient 
processes in the cavity but not fast enough to decouple the mirrors. In this regime the transients caused 
by the mirror motion never die and the cavity never reaches equilibrium. In such a regime the motion 
of one mirror is constantly affected by the other mirror and the delays play crucial role. Dynamic 
processses in LIGO cavities with frequencies between 100 Hz and 10 kHz belong to this regime. 

The third regime takes place when the mirrors move so slow that the laser transients can be ne- 
glected. In this case we can think of the laser circulating in the cavity as being in equilibrium with 
the mirrors. This is the adiabatic regime. In this regime the delays are not important. The adiabatic 
regime takes place in LIGO cavities at frequencies less than 100 Hz. 
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2.7 Self-consistent Field in Adiabatic Regime 

In the adiabatic regime the mirror motion is so slow that we can neglect the delay time in the equation 
for the self-consistent field eq. @ and obtain a solution: 

m = i-w-«**w (23) 

The value of the field in the adiabatic regime is completely defined by the current value of the cavity 
length, z(t). The power in the cavity, which we define as the momentum carried by the forward- 
propagating wave (P = \E\ 2 ) is also a function of the cavity length. The maximum power, P max , 
depends on the mirror parameters and defines the cavity gain: 

P t 2 

^ _ J- max _ ^ a (r . 

mV ~ Pin ~ (l-r a r b r (24j 
The power as a function of length is described by the Airy formula 

p 

p _ 1 max /oc\ 

l + Fsin 2 /^)' 

also known as Airy intensity profile (see [|T7|]). The constant, F, is the coefficient of finesse 

(1 - r a r b ) 2 ' 

The relation between the parameter, F, and the finesse of Fabry-Perot cavity is given in Appendix ^. 
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3 Normal Coordinates 



3.1 Forces due to Radiation Pressure 

The suspended mirrors are interacting with each other through the laser in the cavity. Therefore, the 
degrees of freedom which correspond to the mirrors are coupled. However, we can formulate the 
dynamics in terms of new degrees of freedom which are largely independent of each other. These 
new degrees of freedom are similar to the normal coordinates usually appeared in theory of coupled 
linear oscillators. Before we formulate the dynamics in terms of the normal coordinates we reduce 
the number of fields entering the general expressions for the radiation pressure on the mirrors, equa- 
tions 



The radiation pressure on the end mirror, eq. ( jl5| ) can be expressed entirely in terms of the self- 
consistent field: 

R b (t) = -(l + r 2 b -t 2 b )\E(t-T)\ 2 . (27) 

c 

The radiation pressure on the front mirror is produced by the laser circulating in the cavity and the 
laser incident on the cavity, eq. (JT3]). A constant flux of the incident laser produces constant radiation 
pressure on the front mirror: 

R m = \l + r 2 a -t 2 a )\E m \ 2 . (28) 

c 

Therefore, we can write the radiation pressure on the front mirror as the sum: 

R a (t) = -R b (t-T) + R in + AR(t). (29) 

In this equation the first term implies that the same force which is pushing the end mirror now will be 
pushing the front mirror in the opposite direction after time T. The second term is proportional to the 
power of the incident laser and accounts for the constant radiation pressure from outside the cavity. 
The last term, AR(t), which can be called the excess radiation pressure, accounts for the asymmetry 
between the mirrors and the interference between the incident laser and the laser circulating in the 
cavity. The excess force can also be expressed entirely in terms of the self-consistent field and the 
cavity length: 

AR(t) = -\E(t - 2T)| 2 - - t a r a r b Re{e~ 2lkz ^ E* n E(t - 2T)}, (30) 
c c 

where a = 1 — t 2 — r 2 {r\ — t 2 ), and the asterisk stands for complex conjugation. 

In the adiabatic regime the excess force can be written explicitly as a function of length: 

AR _ _i p a + ^r a r b (r a r b - cos 2k z) 

1 + t Sill KZ 

Numerical estimates with LIGO parameters show that the excess radiation pressure, AR, is roughly 
four orders of magnitude less than the radiation pressure, R b . 



The constant radiation pressure, eq. (|28[), has no effect on the dynamics, but the excess force, 



eq. (P0]l, makes the radiation pressure on the mirrors unbalanced and affects the dynamics. 
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3.2 Normal Modes 



The normal coordinates are the position of the center of mass and the cavity length: 

z( r, = m a x a (t) + m b x b {t - T) ^ 
m a + m b 

z(t) = x b (t-T)-x a (t). (33) 

Note that both normal coordinates are affected by the delay due to the finite time of laser propagation 
in the cavity. 

Similarly, we introduce the control parameters for the new coordinates: 

TT( , m a u a (t) + m b u b (t - T) R in 

m a + m b c(m a + m b )uj$ 

u(t) = u b {t-T)-u a (t)--^. (35) 

To make the subsequent formulas simpler we hid the constant radiation pressure into the new control 
parameters. Also we define the reduced mass in a usual way: 

H = ■ . (36) 

m a + m b 

In normal coordinates the dynamical equations are the equation for the center of mass and the 
equation for the cavity length: 

Z + 1 Z + ulZ = up + AR , (37) 

m a + m b 

1 1 Ai? 

z + 72; + ojn,z = UqU H . (38) 

fj, m a 

We see that the radiation pressure mostly affects one degree of freedom: the cavity length, z(t). This 
is because the pressure on the front mirror is roughly equal and opposite to the pressure on the end 
mirror if we neglect the constant pressure from outside. The motion of the center of mass is affected 
by the radiation pressure only through the excess force. 



3.3 Center of Mass and Observability 

Dynamics of the center of mass of free cavity is defined exclusively by the excess radiation pressure, 
which, in turn, is defined by the cavity length, z(t). As a result, changes in the cavity length excite 
motion of the center of mass. On the contrary, arbitrary motion of center of mass has no effect on the 
cavity length. Therefore we can ignore the dynamics of the center of mass and focus on the cavity 
length. However, few important points must be taken into account before we abandon the motion of 
the center of mass. 

Although small motion of the center of mass (~ 1/xm) has no effect on the detector performance, 
large motion (~ 0.1 mm) will cause variations of the actuator gains. Therefore, the motion of the 
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cavity center of mass must be suppressed to the extent defined by the actuator's range. In particular, 
we must avoid exciting the motion of the center of mass by a control system. In gravitational wave 
detectors the control system suppresses motion of the mirror under seismic disturbances by applying 
forces to the mirrors through the actuators. Therefore the control system can excite the motion of the 
center of mass, unless the control signals are synchronized: 

6u a {t) = -—6u b (t-T). (39) 
m a 

Unlike the cavity length the center of mass is not observable in a strict relativistic sense. To 
obtain information about the cavity's center of mass requires setting up an inertial frame, which is 
impossible for the ground based gravitational wave detectors. This poses a problem of observability 
for control of the Fabry-Perot cavity, which can be stated as follows. The cavity length is known with 
high precision (better than lCT 13 m) but the center of mass can only be measured approximately (up 
to 1/xm). Therefore the cavity length can only be controlled approximately. Since the two degrees of 
freedom are largely independent poor observability and control of the center of mass do not affect the 
dynamics of the cavity length. Therefore, as long as the motion of the center of mass remains small, 
it can be ignored. 



3.4 Cavity Length 

The degree of freedom most important for understanding the response of Fabry-Perot cavity as a part 
of the gravitational wave detector is the cavity length z(t). 

In the following discussion of the cavity length dynamics we neglect the excess radiation pressure, 
AR. To simplify following formulas we make the approximation: 1 + r\ — tf ~ 2 in the eq. ^7] and 
write the radiation pressure on the end mirror as 

R b (t) = -\E(t-T)\ 2 . (40) 

c 

This is a good approximation for high reflective mirrors, such as the end mirrors of LIGO cavities. 

The dynamics of the cavity length is independent of the dynamics of the center of mass and can 
be studied separately. The equations for the cavity length dynamics include only the cavity length and 
the self-consistent field: 

z + 1 z + u 2 z = u 2 u + ^-\E(t-2T)\ 2 , (41) 
E (t) = t a E in + r a r b e~ 2ikz ^ E(t- 2T) , (42) 



The iteration equation ( |42| ) allows us to calculate the amplitude E(t) at discrete times. The time step 
is the cavity round-trip time, 2T. The value of the amplitude between the steps can be found by 
delaying its value at the last iteration. 

These equations are equivalent to the equations of the dynamics of the Fabry-Perot cavity with 
front mirror fixed and the end mirror suspended from wires, provided the mass of the end mirror is 
equal to the reduced mass. 
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4 Equilibrium States 
4.1 Static Solutions 

Equilibrium states of the cavity are static solutions of the equations ( pH] ) and ([4"2|). For static states 
these equations simplify and can be written as 

2 

IMjliz - u) = - \E\ 2 , (43) 
c 

E = i taEm ^ kz - (44) 
1 — r a rf,e 



The complex field, £7, can be eliminated from the equations Q43D and (]44]) and the equilibrium condi- 
tion can be written entirely in terms of the cavity length: 

2/ \ Rmax /ac\ 

Hu {z -u)= 2 ■ (45) 

1 + b sin kz 

Here R max is the maximum radiation pressure on the end mirror which corresponds to the maximum 
power in the cavity: 

_ 2 

Rmax Pmax- (46) 

C 



The equation ( |4"5| ) shows that in equilibrium the radiation pressure is balanced by the pendulum restor- 
ing force. 

Numerical values of the power and the radiation pressure in LIGO Fabry-Perot cavities are given 
in the Table 0} 

Table 1 : Power and radiation pressure in LIGO Fabry-Perot cavities 





P m (W) 


PraaxiW) 


F m ax(N) 


40m LIGO Prototype 


0.25 


160 


1.1 x 10" 6 


4km LIGO 


90 


12000 


7.8 x 10~ 5 



4.2 Strength of Radiation Pressure 

To characterize strength of the radiation pressure we introduce a dimensionless parameter 

a = (47) 

where k = 2n/\is the wave number, ji is the reduced mass and co is the pendulum frequency. The 
parameter a shows how strong the radiation pressure is compared to the pendulum restoring force. 
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This parameter allows us to characterize nonlinearities induced by the radiation pressure quantita- 
tively. As a measure of nonlinearities, the parameter a is similar to the Reynolds number in fluid 
dynamics. 

The parameter a also allows us to establish equivalence of the equilibrium states of different 
Fabry-Perot cavities. Namely, the equilibrium states of a Fabry-Perot cavity with high circulating 
power and heavy mirrors are equivalent to those of a cavity with low power and light mirrors if the 
two cavities have the same a. 

The parameter a has a simple geometric interpretation. From the equilibrium condition, eq. (^3|), 
we see that the maximum increase in the cavity length due to radiation pressure is a/k. Therefore, 
a/n is the maximum number of fringes the mirrors can be pushed apart by the radiation pressure. 

We can think of a as a function of the incident power 

a = ^%P m , (48) 

where G cav is the cavity gain, eq. ([24]). Therefore, the simplest way to change the parameter a is to 
vary the input power of the Fabry-Perot cavity. 

4.3 Equilibrium Condition 

The equilibrium states depend on properties of Airy function: 

1 + t sin 

A discussion of Airy function can be found in ref. [Jl7j]. Some of the properties of Airy function are 
given in the Appendix 

Let us introduce the phase variables for the cavity length and the control parameter: 

= kz, (50) 
rj = ku. (51) 



In terms of the phase variables the condition for equilibrium, eq. (ft5|), becomes 

-(4>-r ] )=A(4>). (52) 
cr 

By solving this equation we can find equilibrium states, 0, for any given value of the control param- 
eter, r\. Even more important is an inverse application of the equilibrium condition. Namely, we can 
turn any state, 0, into an equilibrium state by adjusting the control parameter, rj, according to 

t] = 0- oA{<f>). (53) 

An operation point for the Fabry-Perot cavity is an equilibrium state, at which we want to maintain 
the cavity. We select the operation point by specifying the phase, O , and by adjusting the control 



parameter according to the eq. (53). 
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In LIGO interferometers the operation point is chosen at the maximum of Airy function (0 O = 
nn), where n is integer. To obtain an equilibrium state at these points we must have the control 
parameter set to the following values 



A 



a 



u=- 2 [n-- } . (54) 

Since all the peaks of Airy function are equivalent the different values of the control parameter, 
which correspond to different values of n, are all equivalent. The least absolute value of the control 
parameter corresponds to n = round(cr /ir). Since the control parameter is a linear combination of the 
actuator force and the displacement of the suspension point, eq.([T8|), this condition can be achieved in 
two ways. We can either apply the forces to the mirrors or shift the suspension points. The larger the 
actuator forces the higher the noise in the actuators. Therefore, it is desirable to operate the cavities 
with as less force on the mirrors as possible. 



4.4 Graphical Construction of Equilibrium States 

The equilibrium states are given by the nonlinear equation, eq. which cannot be solved analyti- 
cally. However, we can analyze the equilibrium states using a graphical construction. The graphical 
method for finding the equilibrium states appears in several papers, for example, in ref. [|TE[]. In our 
approach we use the graphical method not only to show approximate locations of the equilibrium 
states but also to derive exact statements about equilibrium and stability of the Fabry-Perot cavity. 

The graphical construction is to plot the functions in the left and the right sides of the equation ( |5"2| ) 
and look for their intersections. The left side of the equation is the pendulum restoring force in 
dimensionless units. It corresponds to a straight line shown in Fig. ^. The line has the slope a -1 and 
the intercept 77. The right side of the equation is the Airy function, which is the radiation pressure 
normalized to unity at its maximum. The points, where the straight line intersects the Airy function, 
are the equilibrium states. In Figure ^| there are three equilibrium states: 0i, 2 and 3 . 



4.5 Condition for Multistability 

A Fabry-Perot cavity is multistable if it has more than one equilibrium state, which belong to the 
same value of the control parameter. To derive a condition for multistability consider a line tangential 
to the Airy profile at some point, O . The slope of the tangent line is A'{<p ), the derivative of Airy 
function at this point. The highest slope of the tangent line corresponds to the maximum value of the 
derivative 

max{A'} « IV3F, (55) 
8 

which is derived in Appendix ^. If the slope of the straight line, which represents the pendulum 
restoring force, is greater than the maximum slope, the line intersects the Airy profile at only one 
point. Otherwise, the straight line and Airy profile can have several common points. Therefore, there 
is a critical value of the parameter a, for which a qualitative change in the equilibrium states occurs. 
This critical value is defined by the maximum slope 

1 8 
maxjA'} 3V3F 
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and depends only on the finesse of the Fabry-Perot cavity. These geometric arguments lead us to the 
condition for multistability. The Fabry-Perot cavity is 

stable, if a < a cr 

multistable, if a > a cr . 

LIGO Fabry-Perot cavities are multistable: a = 3.9, which is greater than the critical value, 
a rr = 0.012. 



4.6 Higher order bifurcations 

At low radiation pressure (a < a cr ) the cavity has only one equilibrium state. If we increase the 
radiation pressure and reach the critical point (er = a cr ) the equilibrium state splits into three states, 
as shown in Fig. |j. Two of the new states are stable and one is unstable. Such a phenomenon is 
common in nonlinear physics and is called the pitchfork bifurcation [|TP|], [|2~0|]. 



If we keep increasing the radiation pressure, at some point another bifurcation occurs and more 
equilibrium states are created. We can continue the process and obtain an infinite number of bifur- 
cations in the Fabry-Perot cavity. The values of the parameter a, at which these bifurcations occur, 
can be found from the graphical construction shown in Fig. [5|. These values are a n ~ irn where 
n — 1, 2, ... is the order of the bifurcation. In terms of power the condition for the nth order bifurca- 
tion to occur is ^ 

Pmax = -C/JUqXtI. (57) 

Exact number of equilibrium states depends not only on the parameter a, but also on the control 
parameter, u. 
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Figure 4: Pitchfork bifurcation 
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Figure 5: Higher order bifurcations 




5 Stability and Hysteresis 

5.1 Multistability of Fabry-Perot cavity 

Fabry-Perot cavity with suspended mirrors is a multistable dynamical system. There are multiple 
equilibrium states in the cavity which correspond to the same value of the control parameter. Some of 
them are stable and some are unstable. A multistable cavity being brought near unstable equilibrium 
state moves to a nearest stable equilibrium state, which belongs to the same value of the control 
parameter. 

Multistability of the Fabry-Perot cavity manifests itself in multi-valued functions describing the 
cavity state. For example, the equilibrium cavity length, z, is a multi-valued function of the control 
parameter, u. This function is defined implicitly: 

z = u + TA(kz), (58) 
k 

and is shown on Fig. ^[ 
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Figure 6: Multistability of cavity length 




A similar multistable response appears if we vary the input power, P in , within a wide range, but 
keep the control parameter fixed. This way we obtain the power in the cavity, P, as a multi-valued 
function of the input power. Such function also is defined implicitly: 



H P 

p ° cav^tn 

" l + Fsm 2 (W^Y 

The plot of this function for a particular value of the control parameter (u = — 2 ^2 ) is shown in 
Fig. [7]. This is a typical multistability curve similar to those, which frequently appear in studies of the 
bistable optical devices [[21]], [|22l]. 

Both plots above have retrograde segments which correspond to unstable equilibrium states of the 
cavity. In general any state on a retrograde part of a multistable curve is unstable [pQ|]. 

The formula which describes the power inside the cavity as a function of the input power, eq. (|5^), 
is equivalent to the corresponding formula for the fixed length cavity filled with a non-linear medium. 
The index of refraction of such medium has a quadratic nonlinearity 

n(E) =n + n 2 \E\ 2 , (60) 

where n and n 2 are the coefficients of the nonlinear index of refraction. A review of the extensive 
research which had been done on the cavities with the quadratic nonlinearity can be found in [ |2T| ] . 
Here we obtain an exact condition for this equivalence. In the equation ( |59| ) the phase is a function of 
the power: 

2k 

<P = ku + M 2 . (61) 

A similar dependence appears in the cavity with fixed length but filled with the nonlinear medium. 
The equivalence is achieved by tuning the parameters according to the equation: 

2k 

n 2 k = o> ( 62 ) 

CjJLUJQ 
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Figure 7: Multistability of power in Fabry-Perot cavity. The parameters are chosen so that the curve 
goes through the LIGO operation point 



incident power, Pin (W) 



where l is the length of the cavity with the nonlinear medium. Therefore, equilibrium states of a 
cavity with suspended mirrors are equivalent to those of a cavity filled with the nonlinear medium if 
the condition, eq. (|62|) is satisfied. This equivalence allows us to extend some results obtained for the 
cavities with the medium to the cavity with suspended mirrors. 

5.2 Hysteresis of Mirror Position 

The instability causes jumps and a hysteresis to appear. In this section we describe the mirror jumps 
and the hysteresis in the mirror position in the Fabry-Perot cavity. 

Assume that one mirror is at rest and the other is moving. Let this motion be so slow that the cavity 
remains in equilibrium. We can produce such motion by a slow change of the control parameter. The 
motion of the mirror is shown on the parameter- state diagram, Fig. [8|. 

Let "A" be a starting point. Assume that we force the mirror to move by changing the control 
parameter from ua to m b . As we increase the value of the control parameter the mirror moves steadily 
from point "A" through point H\ to point H 2 . As the mirror reaches point H 2 a sudden transition to 
point H 3 occurs. This jump happens because both points H 2 and H 3 belong to the same value of the 
control parameter but point H 2 is unstable. After the mirror reaches point H 3 , it resumes its steady 
motion and reaches the point "£>". The sequence of states is 

A^ Fi ^ H 2 ^ B. (63) 

The jump H 2 — > H 3 is accompanied by a rapid build-up of power in the cavity, and, therefore, takes 
time of the order of the storage time to complete. 

In the reverse process, we change the control parameter from ub to ua and thereby force the mirror 
to move in the opposite direction. On the way back the mirror does not retract its path. Instead, it 
follows a sequence 

B -> H 4 ^ Ht^ A. (64) 
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Figure 8: Hysteresis and instability region 




The reverse motion of the mirror also has a jump but at a different place: point H 4 . The transition 
H± — > Hi is accompanied by a rapid drop of power in the cavity. Thus a slow sweep of the mirror 
across the resonance is an irreversible process, known as hysteresis. The complete cycle: 

fli -> ff 2 -> # 3 -> #4 -> #1, (65) 

is the hysteresis loop, which is typical for nonlinear systems. 

A similar hysteresis appears if we keep the control parameter fixed but vary the input power. Such 
hysteresis was observed by Dorsel and collaborators and is described in ref.[Q]. 



5.3 Boundaries of Instability 

The states on the retrograde segments of the multistability curve on Fig. ^| and Fig. |7| are unstable. 
These states form an instability region which is confined by the two tangent lines with the slope a -1 , 
as shown in Fig. |[ The boundaries of the instability region are defined by the equation 

A\<f>) = -. (66) 
a 

There are two solutions of this equation, which are left and right boundaries of the instability region. 
In the large-F approximation the instability region is 

t)'<*-<-2Sp- (67) 

The left boundary corresponds to low power inside the Fabry-Perot cavity. The right boundary 
is located very close to the peak of the Airy function and corresponds to high power. For example, 
in the 4km LIGO interferometers the right boundary of the instability is 0o = — 7.5/irad. Thus the 
separation between the right boundary and the peak is very small, roughly a thousandth of the width 
of Airy function. 
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6 Hamiltonian Dynamics 



6.1 Conditions for Adiabatic Motion 

Motion of the mirrors belongs to the adiabatic regime only if it is slow; the characteristic time of 
such motion must be much greater than the cavity storage time. We define the characteristic time 
differently for small- and large- amplitude mirror motions. 

We consider the amplitude small if it is less than the width of the resonance. Usually such motion 
is oscillatory. Then the characteristic time is equal to the period of these oscillations, T osc . Therefore, 
the motion is adiabatic if T osc 3> r. 

Large-amplitude mirror motion takes place when the changes in the cavity length are greater 
than the width of resonance, x w . During such motion the relative velocity of the mirrors, v, can be 
considered constant within the time the mirrors pass through a resonance. Then the characteristic 
time is the time it takes for the length to change by the width of one resonance: x w /v. Therefore, the 
motion is adiabatic if x w /v 3> r. This condition can also be seen as a requirement on the relative 
velocity of the mirrors. Namely, the large-amplitude mirror motion is adiabatic if v <C v a d, where v a d 
is the adiabatic threshold velocity: 

v ad = — ^ — — . (68) 

t nr Lq 

The threshold velocity depends only on the parameters of the cavity and the wavelength of the laser. 
For LIGO Fabry-Perot cavities v a d = 7.4 x 1CT 7 m/s. 



6.2 Effective Potential 

In the adiabatic regime the field inside the cavity, E(t), is completely defined by the length, z(t), see 
eq. ([23]). Therefore, we can eliminate the field and consider only the dynamics of the length. The 
equation for the length dynamics becomes 

z + jz = -uJq(z - u) + -R max A(kz), (69) 



where fi is the reduced mass and -R max is the maximum radiation pressure defined by eq. (|4q). This 
equation describes the relative motion of the mirrors driven by the pendulum restoring force and the 
radiation pressure. A similar equation for the Fabry-Perot cavity with one suspended mirror appeared 
earlier in refs. [§, QTTJ, pTTp . 

The pendulum restoring force is generated by the usual quadratic potential 

V p {x) = ^ul{z-u)\ (70) 

which depends on the control parameter, u. The force due to the radiation pressure is generated by a 
potential 

arctan[ v / F + ltan(fcz)] 

ky/F + 1 

which is derived in the Appendix 0. 
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The sum of the two forces is generated by the effective potential 

V eS (z) = V p (z)+V iad (z). (72) 

Thus the relative motion of the mirrors interacting with the electro-magnetic field in the cavity is 
equivalent to the motion of a particle with the reduced mass in the field of the effective potential. 
Such representation is valid only in the adiabatic regime, when the delay time can be neglected. 
Beyond the adiabatic regime the delay time becomes important and the field depends on history of 
mirror's motion within the storage time. 

A remarkable property of the adiabatic regime is that the dynamics is hamiltonian. In other words, 
we can define a total energy for the motion in the effective potential: 

S = ^ + V cS (z). (73) 

Thus the Fabry-Perot cavity in the adiabatic regime is a hamiltonian system. Since there is friction 
(7 > 0) the energy is not conserved: 

£ = _ 7/i i 2 < 0. (74) 

However, the friction in suspension systems of the gravitational wave detectors is very small. Con- 
sequently the energy changes very little (less than a part per million) over a period of the pendulum. 
Therefore, on time scales of the order of few seconds and even minutes we can think that the energy 
is conserved. 

Note that if there is no friction (7 = 0) the energy of the motion in the effective potential, eq. (f73|). 
is conserved but the combined energy of the mirrors and the electro-magnetic fields in the cavity may 
not be conserved. The cavity, being multistable, may release or absorb energy as it moves from an 
unstable to a stable state. 



6.3 Global Properties of Effective Potential 

The effective potential as a function of cavity length depends essentially on three parameters: F, a 
and r\. This can be seen if we write the effective potential in terms of dimensionless parameters: 



eff 



yuJl \ ((f) — i]) 2 arctan( y/F + 1 tan ( 
a — 



A' 2 



VF + 1 



(75) 



The overall constant provides correct dimensions to the potential, but otherwise is irrelevant. The 
parameter a and the control parameter r] can be dynamically adjusted, but the coefficient of finesse, 
F, cannot be. Therefore, we only consider how the effective potential changes with respect to arbitrary 
changes in the two parameters: a and 77. 

Global properties of the effective potential are defined by one parameter only, a. For a < a cr the 
potential has only one minimum and for a > a cr there are several local minima. Effective potential 
for two different values of o is shown in Fig. [5]. The plots correspond to stable and multistable Fabry- 
Perot cavities. 

The effective potential as a function of two parameters establishes a connection between the dy- 
namics of the Fabry-Perot cavity in the adiabatic regime and the dynamical systems of the catastrophe 
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theory [|Z3|]. The catastrophe theory studies singularities (abrupt changes) in equilibrium states of a 
dynamical system caused by a continuous change of its parameters. The pitchfork bifurcation, as 
described above, is one of the examples of such singularities. 



6.4 Minima and Maxima of the Effective Potential 

The control parameter, r\, allows us to fine tune the effective potential. By changing the control 
parameter we can shift the positions of local minima and maxima of the potential. To illustrate this 
consider the following example. Let the operation point be 0i as shown in Figure [K]. 

To make this state an equilibrium we have to adjust the control parameter according the eq. (|51|). 
This is equivalent to solving the following equation: 

dV eS 



0. 



Since at the point, 



the effective potential has a local minimum: 

d 2 V eS 



> 



(76) 



(77) 



this state is stable. From the figure we see that if we choose the operation point at <p 2 , we obtain 
the same control parameter and arrive at the same effective potential. However, at the point <p 2 the 
effective potential has a local maximum: 



d 2 V cS 



< 0. 



(78) 



4>=4>2 



Therefore, the equilibrium state at 2 > is unstable. 

Thus, stable equilibrium positions correspond to local minima and unstable equilibrium positions 
correspond to local maxima of the effective potential. By adjusting the control parameter we can turn 



any state into a stable equilibrium if that state is outside the instability region, eq.(p7|). 
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7 Resonance of Cavity Length 
7.1 Resonant Frequency 

A stable equilibrium state of the cavity corresponds to a local minimum of the effective potential. Let 
the minimum of the effective potential be at z . The small-amplitude mirror motion near this equi- 
librium state produces harmonic oscillations of the cavity length. The frequency of such oscillations, 
fl , is defined by the curvature of the potential at the minimum: 



(79) 



The oscillations of length due to motion of a single mirror, confined in one of the minima of the 
effective potential were observed by Dorsel et al [^|, [7|]. 

Since there are multiple equilibrium states, there are multiple frequencies for the oscillations of 
the cavity length. Each minimum of the effective potential has its own frequency. These are the 
resonant frequencies of the cavity length, if the cavity is driven by the external sinusoidal force. At 
high radiation pressure these frequencies substantially deviate from the pendulum frequency, c^o- We 
find that the frequencies can be expressed entirely in terms of the parameter a and the location of the 
minimum: 

^(0 o )=^[i_ a A'(0 o )], (80) 

where O = kz and A'((j) ) is the derivative of the Airy function (Appendix 0). The plot of f2 /27r 
as a function of the phase O with LIGO parameters is shown in Fig. [TT]. 
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Figure 11: Frequency of length oscillations as function of operation point 
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The frequency reaches maximum at the right inflection point of the Airy profile. The maximum 
frequency is 

Winfl) = ouo (l + — V . (81) 



In LIGO Fabry-Perot cavity the maximum frequency is 15.7 Hz. 



7.2 Adjustment of Resonance Frequencies 

The resonances of cavity length make the Fabry-Perot cavity a narrow-band transducer between me- 
chanical and optical signals. This observation was made by Dorsel and collaborators They also 
suggested that a resonant frequency can be tuned by changing the power of the incident laser. The 
adjustment of the incident power, which are equivalent to the adjustments of the parameter a is not 
the only way to tune the resonant frequency. We can also vary the control parameter. By changing the 
control parameter we can shift the minimum of the effective potential and thus change the resonant 
frequency corresponding to that minimum. 

On the other hand we may require that the operation point remains fixed as we tune the resonant 
frequency. This can be done if we vary the two parameters, a and 77, simultaneously so that the 
equilibrium condition is preserved. For this we need to accompany any change 5a by the change in 
the control parameter as follows 

8r) = -A((f) ) 5a. (82) 

This way we change the curvature of the effective potential at the minimum without changing the 
location of that minimum. 

The result is the resonant frequency as a function of the control parameter 



Oo - V) 



(83) 
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Figure 12: Fine tuning of effective potential 




The fine tuning of the effective potential and the corresponding frequencies are shown in Fig. [T2 



7.3 Plateaus of the Effective Potential 



By tuning the parameters, o and rj, we can make the effective potential more and more flat at one of 
its minima and even form a plateau: 

' 2T/ eff 



d 2 V P , 



dz 2 



0. 



(84) 



One way of creating the plateau is to tune the control parameter, r\. In this method we flatten the 
potential by bringing the minimum closer to the nearest local maximum. This method of obtaining 
the plateau by merging the local minimum with the nearest local maximum has a drawback. Namely, 
in this method we bring the operation point closer and closer to the instability region, eq. (]57|). 

Another way of forming the plateau at the minimum is to start with lower power (a < o cr ) and 
gradually increase the power until the bifurcation occurs (a = o cr ). In this method we begin with a 
stable cavity: the effective potential has only one minimum. We have to adjust the control parameter 
so that the operation point coincides with this minimum: 



d 3 V e 



off 



dz 3 



0. 



(85) 



By increasing the incident power o — > a CT , we make the potential more and more flat at the minimum. 
At the bifurcation point (a = a cr ) the plateau is formed. 
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7.4 Period of Anharmonic Oscillations 



The oscillations of length in almost flat potential are anharmonic and can have very large period. At 
the critical point the potential is of fourth order: 

V eS (z) = V cS (z ) + ^V}*\z )(z - z )\ (86) 
The period of the anharmonic oscillations in the fourth order potential is 



6/i 



*v}*\z )\ a 



i 



2 1 

(87) 



where a is the amplitude of the oscillation. Since the amplitude of mirror motion and correspondingly 
the amplitude of the length oscillations is very small (of the order of a few /im), the period can be 
very large. 

In practice it is unlikely that the cavity is maintained exactly at the bifurcation point. More likely 
the long period oscillations are obtained slightly below the critical point a < a cr . In this case the 
potential is not quartic, nonetheless the oscillations have long period, which can be found from a 
general expression 

dz 



Tosc = p» / - -t , (88) 

where z\ and z<2 are the turning points of the oscillations. Both equations ( ]8"7| ) and (|8"%|), appear in 
books on Mechanics, for example in ref. [|24|]. 

In Table [Z] we show numerical values of frequencies of the oscillations in the nearly flat potential 
for LIGO Fabry-Perot cavities. The flatness of the potential is characterized by the deviation of the 
parameter a from its critical value. The table shows that we can achieve period substantially greater 



Table 2: Period of anharmonic oscillations in almost flat potential 



amplitude (m) period (s) 

a a = 0.96 a cr a = 0.98 o cr a = a cr 

HP L34 L34 L33 - 

10~ 9 2.31 2.36 2.41 

10~ 10 6.39 8.74 23.2 

10" 11 6.68 9.44 234.2 



than the pendulum period (2tt/ljq w 1.3 s). This method of obtaining very long period by taking 
advantage of the radiation pressure is an interesting alternative to the very long suspension. 
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8 Implications for Control System 



8.1 Pound-Drever Locking Servo 

The operation point of the cavity is a stable equilibrium length which corresponds to a minimum of 
the effective potential. Due to the ambient seismic motion the actual state of the cavity is constantly 
moving near the operation point. Often the minimum of the effective potential is not deep enough to 
confine the state. In such case a negative feedback control system (servo) is used to keep the cavity 
state within the minimum. The dynamics of the cavity with the control system is described by the 
same equations as before but the control parameter, u, becomes a function of the cavity length. In 
LIGO interferometers the control function, u(z), is proportional to the Pound-Drever signal, eq. (fTl]): 

u{z) oc lm{E(t)}. (89) 

Such control function depends on the cavity length z and does not depend on the cavity center of 
mass. 

There is no explicit analytic expression for the function u(z), given by eq. (|8"9]), which is valid for 
all dynamical regimes. However, in the adiabatic regime we can obtain a simple formula 

, , / G \ sin 2kz 

"W = -UJ i + + (90) 

where G is the gain and u is the dc-bias of the control system. The minus sign in eq. ( |9T| ) accounts 
for the negative feedback. 

If the deviation of the cavity length from the nearest fringe is much less than the width of the 
resonance the control signal becomes linear: 

u(z) « -G^z-n^j +u , (91) 

where n is the order of the fringe. Since the amplitude of motion of the mirrors is small (of the order 
of a micron) the fringe order number is a small number: n = 0, ±1, ±2, .... 



8.2 Multistability in Presence of Servo 

Since the control law eq. ( |9T!| ) is a periodic function of the cavity length the cavity with the servo is 
multistable. To study multistability in dynamics of the cavity with the control system we introduce 
the static servo potential 



V svo (z) = -/jwq J u(z r )dz' 
o 

oc ln(l + Fsin 2 kz). 



(92) 



(93) 



Therefore, in the adiabatic regime the dynamics of the cavity with the servo can be described by the 
combined potential 



k 2 



— a 



ar ct an ( a/F+T t an < 



G_ 
2F 



In AO) 



(94) 
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where the dc-bias is conveniently assigned to the pendulum potential (tj = ku ). 

For large- amplitude mirror motion the control function, u(z), is nonlinear and, therefore, the servo 
acts as a nonlinear device. The effect of the servo depends on the gain, G, and can be made stronger 
than the effect of the radiation pressure. Fig. O shows the sum of the radiation pressure potential and 
the servo potential for three different values of the servo gain: G = 10 2 , 10 3 , 10 4 . On the first plot 



Figure 13: Combined radiation pressure and servo potential 
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the sum of the two potentials looks like a staircase which is characteristic of the radiation pressure 
potential. Here the effect of the servo is negligible compared to the effect of the radiation pressure. 
On the second plot the "staircase" shows dips created by the servo. Here the effect of the servo is 
comparable to the effect of radiation pressure. On the third plot the dips are bigger than the steps, and 
the servo dominates the radiation pressure. 

8.3 Equilibrium and DC-bias 

Stable equilibrium states of the cavity with the servo are the minima of the combined potential: 

^-(V eS + V sm ) = 0. (95) 
az 

This condition can also be seen as a requirement on the dc-bias. Namely, a cavity state, z , becomes 
an equilibrium if the dc-bias is set according to the equation: 

uo = z - ^ ^1 - — sin2kz S j A(kz ). (96) 

From this equation we see that, in general, the dc-bias is connected to the gain of the servo. There- 
fore, changes of the servo gain should be accompanied by the corresponding changes of the dc-bias. 
Otherwise the equilibrium condition will be lost and the equilibrium state will change. There is one 
exception to this rule. In the case when the operation point is chosen at the peak of the Airy function 
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(zo — the dc-bias becomes independent of the gain: 



A a 



(97) 



Un = n — . 

2 k 



Note the dependence of the dc-bias on the fringe number, n. 

8.4 Stability Provided by Servo 

The negative feedback control system can turn an unstable state of the cavity into a stable one. It can 
also provide greater stability margins to the already stable cavity. These effects of the servo can be 
analyzed with the help of the servo potential. 

The frequency of the cavity length oscillations in the presence of the servo is 



By adjusting the gain we can vary this frequency in a wide range. We can also turn the unstable state 
(fig < 0) into a stable one (Qq > 0). In particular we can make stable any state of the cavity within 
a width of the fringe (|0 O | < <t>w)- This can be done by choosing a sufficiently high servo gain. For 
LIGO parameters such gain is G « 1800. 




"^(Kff + ^VO) 

o 2 {l - aA'ifo) + G[l - (F + 2) sin 2 o ]^ 2 (0o)} , 



(98) 



(99) 
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9 Conclusions 



The suspended mirrors in the Fabry-Perot cavity are described as a pair of harmonic oscillators driven 
by nonlinear forces of the radiation pressure. The oscillators are not independent; they interact with 
each other through the laser circulating in the cavity. The normal coordinates for the coupled oscil- 
lators are the position of their center of mass and the cavity length, both affected by the delay in the 
cavity. The radiation pressure has almost no effect on the motion of the center of mass, but affects 
strongly the dynamics of the cavity length. To understand universal properties of the Fabry-Perot 
cavity as a nonlinear dynamical system we introduced two parameters. These parameters determine 
qualitative behavior of the dynamics of the cavity length common to all Fabry-Perot cavities with 
suspended mirrors. The first parameter determines whether the cavity is multistable or not. There is a 
critical value for this parameter, at which the bifurcation occurs, and the cavity becomes multistable. 
The second parameter is the generalized control force. It determines the locations of the equilibrium 
states. If the control parameter varies within a wide range the equilibrium state of the cavity follows a 
hysteresis loop. We analyze stability of the equilibrium states and identify the instability region. The 
instability is explained in terms of the effective potential: stable states correspond to local minima of 
the effective potential and unstable states correspond to local maxima. Each minimum of the effective 
potential defines the resonance frequencies of the cavity length oscillations. These resonances make 
the cavity a narrow band detector of any disturbances which affect the optical path of the laser in 
the cavity, such as mechanical vibrations of the mirrors. The resonance frequencies can be tuned by 
changing both of the parameters. 

The nonlinear dynamics, described in this paper, provides grounds for building a control system 
for the Fabry-Perot cavities in the interferometric gravitational detectors. Our results can also be used 
in the studies of the Fabry-Perot cavities near the quantum limit. 
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A Parameters of LIGO Fabry-Perot Cavities 

The following table shows the transmissivities and the losses for the 40m and 4 km LIGO mirrors. 



Table 3: Transmissivity and reflectivity of LIGO mirrors 



Mirror # 


parameter 


40m Prototype 


4km LIGO 


input mirror 


t a 


0.0758 


0.1732 




r a 


0.9971 


0.9849 


end mirror 


h 


0.0035 


0.0045 




n 


0.999944 


0.999965 



Other parameters for LIGO cavities are listed in the following table 



Table 4: Parameters of 40m and 4km LIGO cavities 



Parameter (units) 


Symbol 


40m Prototype 


4km LIGO 


laser wavelength (m) 


A 


5.14 x 10~ 7 


1.06 x 10~ l 


cavity length (m) 


L 


38.5 


4000 


delay time (s) 


T 


1.28 x 10~ 7 


1.33 x 10- 


storage time (s) 


T 


8.60 x 10~ 5 


1.74 x 10- 


free spectral range (Hz) 


/fsr 


3.89 x 10 6 


3.75 x 10 4 


cavity low-pass freq. (Hz) 


fcav 


1.85 x 10 3 


91.2 


cavity gain 


G cav 


645 


130 


coefficient of finesse 


F 


4.47 x 10 5 


1.71 x 10 4 


phase width (rad) 


4> w 


1.49 x 10- 3 


7.64 x 10- 


displacement width (m) 


x w 


1.22 x 10~ 10 


1.29 x 10-' 


Finesse 


Finesse 


1051 


205 


mirror mass (kg) 


m a 


1.53 


10.8 


pendulum frequency (Hz) 


uj /2ti 


1.0 


0.77 
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B Derivation of Field Equations 

We approximate the laser by a plane electro-magnetic wave with electric field, f(x,t), which is a 
solution of the 1 -dimensional wave equation: 

1 d 2 f d 2 f 

= 0. (100) 

c 2 dt 2 dx 2 

Let fi(x, t) be the electric field corresponding to the right-moving wave in the cavity and f 2 (x, t) be 
the electric field corresponding to the left-moving wave. We define the complex amplitudes at the 
reference planes, located at x = and x = L . The amplitudes of the right-moving wave, E\ and E[, 
are defined by the following equations: 

A(0,t) = E 1 (t)e^ t , (101) 
A(Lo,t) = E[{t)e^. (102) 

Similarly, we define the amplitudes, E 2 and E' 2 , for the left-moving wave: 

f 2 (0,t) = E 2 (t)e^, (103) 
f 2 {L ,t) = E 2 (t)e wt . (104) 

Since the wave functions satisfy the wave equation their values at different locations can be obtained 
by delay: 

/i(Lo,t) = A(0,t-L /c) 7 (105) 
f 2 (0,t) = f 2 (L ,t-L /c). (106) 

Therefore, the amplitudes at the two reference planes are related to each other: 

E[{t) = E x (t-Ulc)e- lkL \ (107) 

E' 2 (t) = E 2 (t- L /c)e- ikLo . (108) 

If we assume that 

L = N^, (109) 



where N is an integer, we obtain the equations (0-0)- 

The equation for the reflection off the end mirror can be obtained from the condition of continuity 
of the wave function at the mirror surface: 

f 2 (L + x b ,t) = -r b fi(L + x b ,t). (110) 

Therefore, at the reference plane x = L : 

f 2 (L ,t + x b /c) = -r 6 /i(L ,t- x b /c). (Ill) 

This equation can now be written in terms of the amplitudes: 

E 2 (t + x b /c) = -r b E[{t - x b /c)e- 2%kx \ (112) 

If we neglect the small time delay due to the mirror motion, x b /c, we obtain the equation Other 
equations can be derived similarly. 
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C Some properties of Airy function 



Here we present few basic properties Airy function, which are used throughout this paper. The Airy 
function, 

A((f>) — - ; ^7~r~2~T ) (H3) 



1 + Fsin^ 

is a postitive, periodic function of <p with the period ir. The Airy function depends on one parameter, 
F, a positive number called the coefficient of finesse. Fig. [13] shows Airy functions with different 
coefficients of finesse. 

Figure 14: Airy Profile 




-4-3-2-1 1 

phi 



The half-width, (f> w , is the phase, at which the Airy function equals to half of its maximum value 

1 1 



arcsm 



(114) 



The full width is equal to 2<p w . Finesse is the ratio of the separation between the peaks to the full 
width of the peak 

Finesse = « ^-Vf. (115) 

A slope of a tangent line to Airy function at point (f> is defined by the derivative of Airy function 

F sin 20 



m) = - 



1 + Fsin 2 (b) 2 ' 



(116) 



The inflection points of Airy function can be found by solving the equation: A" ((f)) = 0. They are 
located at ±0 in f, where 



cos20 inf = — (V9F 2 + AF + 4 - F - 2) . 



(117) 
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The greatest slope the tangent line achieves at the left inflection point </> = —4>mi- This maximum 
slope is 

1/2 



\{F + 2V9F 2 + AF + 4 - §F 2 - 2F - 2 

A'(-<f>mf) = ~ ; — — ri ; 

(i V9F 2 + AF + 4 - \F - ' r 



(118) 



The expression for the inflection points, eq. fll 17D , and the maximum slope, eq. ( ]1 18D , are quite 



complicated, therefore, it useful to find simple approximate formulas for these quantities. For high 
finesse cavities (large F- limit) such approximations can be easily found: 

' (119) 



3F 

A'i-M w lV3F. (120) 

o 

Any straight line, which has a slope greater than the maximum slope, intersects the Airy function at 
only one point. If the slope of the line is less than the maximum slope the line can intersect the Airy 
function at more than one point. 

The integral of Airy function can be found in terms of elementary functions 

A{#)0 = arCtan( jJ^ tan0) - (12D 

u 

This expression is used in defining the potential for the radiation pressure in the adiabatic approxima- 
tion. 
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